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ABSTRACT 

Context. Measurements of the magnetic field in the interplanetary medium, of the sunspots area, and of the heliospheric current sheet 
position, reveal a possible North-South asymmetry in the magnetic field of the Sun. This asymmetry could cause the bending of the 
heliospheric current sheet of the order of 5-10 deg in the southward direction, and it appears to be a recurrent characteristic of the 
Sun during the minima of solar activity. 

Aims. We study the North-South asymmetry as inferred from measurements of the deflection of polar coronal hole jets when they 
propagate throughout the corona. 

Methods. Since the corona is an environment where the magnetic pressure is greater than the kinetic pressure {[} <K 1), we can 
assume that magnetic field controls the dynamics of plasma. On average, jets during their propagation follow the magnetic field lines, 
highlighting its local direction. The average jet deflection is studied both in the plane perpendicular to the line of sight, and, for a 
reduced number of jets, in three dimensional space. The observed jet deflection is studied in terms of an axisymmetric magnetic field 

model comprising dipole {g i), quadrupole {gi), and esapole (gA moments._ 

Results. We measured the position angles at 1 R Q and at 2 R e of the 79 jets from the catalogue of lNistico et all d2009h . based on the 
STEREO ultraviolet and white-light coronagraph observations during the solar minimum period March 2007-April 2008. We found 
that the propagation is not radial, in agreement with the deflection due to magnetic field lines. Moreover, the amount of the deflection 
is different between jets over the north and those from the south pole. Comparison of jet deflections and field line tracing shows 
that a ratio gilgi — -0.5 for the quadrupole and a ratio gi/gi — 1.6 - 2.0 for the esapole can describe the field. The presence of a 
non-negligible quadrupole moment confirms the North-South asymmetry of the solar magnetic field for the considered period. 
Conclusions. We find that the magnetic deflection of jets is larger in the North than in the South of the order of 25-40%, with an 
asymmetry which is consistent with a southward deflection of the heliospheric current sheet of the order of 10 deg, consistent with 
that inferred from other, independent, datasets and instruments. 
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1. Introduction 


The solar corona is an environment highly structured by the 
strength and the topology of magnetic fields. Even during so¬ 
lar minima, the corona is far from being a quiet region but it 
evolves on many time scales, including the solar cycle. The ob¬ 
servations of the corona in EUV and X-ray wavelengths reveal 
bright (dense) regions, coinciding with the presence of active re¬ 
gions, and extended dark (void) areas, named coronal holes. In 
the same way, white-light observations obtained from corona- 
graphs show the presence of ray-like features at higher latitudes, 
and helmet streamers at middle and equatorial latitudes. In as¬ 
trophysics, knowledge of the solar magnetic field comes mainly 
from measurements of splitting atomic lines due to the Zeeman 
effect. This method can be applied to the radiation coming from 
the photosphere, allowing to estimate the vector magnetic field, 
but there is not much possibility to have direct measurements of 
the magnetic field of the corona. Indirect estimates of the coro¬ 
nal magnetic field are obtained from extrapolations techniques 


(PFSS, NLFF) dWiegelmann & Sakuraill2012l). bv radiophysi- 
cal_methods and coronal seismology (e.g. lNakariakov & Ofmanl 
1200 ll) . 


There are several observations suggesting a North-South (N- 
S) asymmetry of the solar magn etic field during solar minima 
(e.g., see the discussion in lErdos & Baloghll2010l) . Early evi¬ 
dences came from direct measurements of the photospheric mag¬ 
netic field by magnetograms, which were extrapolated to the so¬ 
lar wind source surface by iHoeksemal (119951) . who showed that 
the magnetic field strength in the Sun’s south polar cap was 60% 
larger than in the north one during some solar minima. Then, the 
Ulysses spacecraft provided new insights of the N-S asymme¬ 
try: thanks to its orbit, nearly perpendicular to the ecliptic plane, 
Ulysses explored high latitudes regions of the heliosphere, and 
during its passage close to the Sun, the so-called fast latitude 
scan at about 1.4 A U, gave us measuremen ts of the interplan¬ 
etary magnetic field dErdos & Baloghll 19981) and particle data. 
Indeed, an indication of N-S asymm etry comes from the global 
distribution of the solar wind speed dTokumaru et al.li2015h . the 
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latitudinal gradients of e nergetic particle fluxes dSimnson et al.1 
Il996h iHeber et~aD[T99 6) which show an unbala nce of fluxes be- 
tweeen north and south of 6% -15%. Moreover. Erd os & Baloghl 
( '201 Oh studied and compared magnetic field data of Ulysses dur¬ 
ing the first latitude scan (coincident with the minimum of solar 
cycle 22 in 1994-95) and the third one (minimum of solar cy¬ 
cle 23 in 2007-08), and shows that the radial component of the 
magnetic field (normalised to 1 AU) of the South hemisphere is 
greater by a factor 1.12-1.21 than that of the North. In addition 
they found a southward shift of the heliospheri c cu rrent sheet 
(HCS ) of the order of 3-5 deg (see also IVirtanen & Mursulal 
l20Toh . A previous limi t of about 7 deg was found by Mursula & 
Hiltula T2004l ) based on heliospheric magnetic field observations 
at 1 AU. As of now, the actual values of the southward shift of 
the heliospheric current sheet remain not well known. 

A possible model for explaining the N-S asymmetry is to 
consider the global magnetic field as the contribution of sev¬ 
eral multipole components, and ascribe the mismatch between 
the north and the south magnetic fields to the contribution of 
the quadrupole moment, as suggested bv Bra vo & Gonzalez- 
Es parza (j200d) an dlMursula & Hiltulal (120041) . Indeed, Fig. 1 
of iBravo & Gonzalez-Esnarzal d200(iirand Bravo & Gonzalez- 
Esparza ( 200llf . which compares the dipole, quadrupole and 
esapole structures, shows how this asymmetry could be gener¬ 
ated: if we look at the direction of magnetic field for the dipole 
and the quadrupole during two consecutive solar minima, we 
note that they are opposite in the north pole, and concordant in 
the south pole. The direction of the magnetic field of the esapole 
is always concordant with that of the dipole at the poles. This 
means that, due to the action of the quadrupole moment, the 
northern magnetic field is weakened, while the southern field is 
enhanced, causing the evident asymmetry and the shift of the 
heliospheric current sheet in the southward direction. 

Observations of non-radial white-light coronal streamers 
(Wan 311996lb or deflection of EUV polar plumes (de Patoul 
et al. l2013bl al) allowed to study the possible configuration of 
the global and polar magnetic field, respectively. On the other 
hand, coronal jets are transient density enhancements of plasma 
that can highlight the local structure of the magnetic field. In¬ 
deed, hot jets at poles observed in EUV or X-rays are naturally 
explained as the result of magnetic reconnection between emerg¬ 
ing flux with the open magnetic field: the tension force of the dis¬ 
connected field lines compres ses the plasma and pushes it away 
along the magnetic field lines dYokovama & Shibatall 1995li . 

The aim of this work is to estimate the configuration of 
the solar magnetic field, and consequently the N-S asymmetry, 
starting from the latitudinal deflection of polar jets observed by 
STEREO spacecraft. Section 2 presents the analysis of the de- 
flectionofthe - £oto_coronaljets i JromJhe_catalogue -i o£Nistic6 
et al. d2009h , which regards jets observed from March 2007 to 
April 2008, which roughly corresponds to the solar minimum, 
between the end of solar cycle 23 and the beginning of 24. The 
analysis is performed in 2D and 3D space. Modeling and ex¬ 
trapolation of the coronal magnetic field structure is explained 
in Section 3, and conclusions are given in Section 4. 


2. Polar jet deflection measurements 

2.1. 2D analysis 

We investigate the influence of the large scale coronal magnetic 
field on the motion of 79 jets through the corona observed with 
STEREO during Mar ch 2007-April 2008, which are catalogued 
in lNistico et al.l ( 2009 1. During this period, the angular separation 


between the two STEREO spacecraft increased from 2 to ~48 
deg. We exploited data from EUVI and COR1 instruments of the 
SECCHI package. Our method consists of calculating the posi¬ 
tion angle (PA), i.e. the angular displacement from the north axis 
to the jet position in counterclockwise direction, at two fixed dis¬ 
tances from the solar centre: 1 R Q , that corresponds to the solar 
limb, and 2 R e (see Fig. [I]). Since the EUVI images cover the full 
disk up to 1.4 R q , they are appropriate for measuring the PA at 
1 Rr.j (8 euvi ), that in most cases corresponds with the position of 
the jet footpoint. COR1 images have a FOV within 1.4-4R 0 , but 
the exact limits vary due to the offset of the occulter with respect 
to the Sun center. In order to clearly identify the jet in COR1, 
we took PAs at 2 R 0 (0 COR1 ), somewhat above the boundary of 
the occulter from the Sun centre. FITS file data were prepped 
within the SolarSoft (SSW) environment by secchi_prep. pro, 
which allows to calibrate and apply corrections to the images, 
including rotation of the axis towards the solar North. Informa¬ 
tion about the Sun’s centre in units of pixels are retrieved for 
EUVI from the keywords CRPIX1 and CRPIX2 of the header, 
for COR1 we used the function wcs_get_pixel. prcQ. the pixel 
size in arcsec is given from CDELT1/2, and the solar radius 
from RSUN. Circumferences defining the distance of 1 and 2 
R q from the centre are over-plotted in the images and the PA 
for each single jet is visually determined and measured in deg 
either from STEREO A or B, according to which spacecraft of¬ 
fers a better view in both instruments. We assume an error of 
A/ =10 pixels in the location of the jet (this apparently large er¬ 
ror can include some other instrumental errors, jitter, etc.) with 
both instruments. This can be easily converted in terms of angu¬ 
lar displacement, by considering the pixel size of each detector 
(A p ix) and the radial distance from the centre, which are not con¬ 
stant through the period of observations. Thus, a typical error in 
EUVI is A d EUVI = Al A pix /r ~ 1 deg, ( with A pix * 1.6" and 
r ~ 950"), and in COR1 Ad C0Rl ~ 2 deg (with A pix ~ 7.5" and 
r ~ 1900"). The measurements of the PAs for the 79 jets from 
COR1 are plotted as a function of PAs from EUVI in Fig. [2]for 
the North pole (left) and the South one (right). The somewhat 



Fig. 2. Plots of the EUVI PAs (horizontal axis) vs the COR1 PAs (ver¬ 
tical axis) for jets occurring at the north (left) and at the nouth polar 
coronal hole (right). The red dashed line represents the bisector of the 
plane and ideally should trace equal PAs between the EUVI and COR1 
FOV (radial propagation); the green dashed line fits the data points. 

Events rooted significantly away from the poles in EUV tend to have an 
even larger deviation in COR1. 

larger number of events in the north coronal hole with respect 
to the south (45 against 34, respectively) is due to the orbital 

1 The keywords CRPIX1/2 for COR1 return the centre of the occulter. 
Coordinates information can be retrieved with the Word Coordinate 
System (WCS) routines in SSW starting from f itshead2wcs. pro. See 
http://hesperia.gsfc.nasa.gov/ssw/gen/idl/wcs/wcs_tutorial.pdf 
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Fig. 1. Image of the Sun at 304 A (left) from EUVI and from COR1 in white light (right). PAs are from the north solar axis in the counterclockwise 
direction (see the reference frame overlapped on the images) at 1 R 0 in the EUVI FOV and at 2 R a in the COR1 FOV. 


features of STEREO, which allow a betterview of the northern 
region in the investigated time period dNistico et al.ll2009l) . In 
addition, during the period of observations the Sun was at mini¬ 
mum. The jet angular positions, represented by circles, are fitted 
by a linear function separately for the north/south pole coronal 
hole (NPCH/SPCH), according to the equation: 

(f ORl = a e EUVI + b (i) 

where a is the slope of the line and b is the constant term of 
the linear fit. We performed linear fits with two different IDL 
routines: the first is a fit unweighted by measurement errors with 
linfit.pro, which performs regression of the data points along 
the vertical distance (dashed green lines in Fig. [2]). and th e sec- 
ond with linfitex.pro of the MPFIT package (iMarkwardtl 
l 2009l) . which performs full-Cartesian regression by taking into 
accounts errors in both x and ij variables (yellow line). Both the 
fits are almost coincident. The red line is a reference line corre¬ 
sponding to the case when jets had the same PAs at 1 and at 2 R Q , 
implying radial propagation. Although some scatter of the data 
points is present, an overall trend is evident: jets having small 
PAs (if in the NPCH) or small displacements from the south po¬ 
lar axis (if in the SPCH) in EUVI FOV show a small deviation 
in the COR1 FOV (they are near or on the red line, representing 
events that have the same PA when seen from both instruments); 
jets having large PAs (if in the NPCH) or large displacement 
from the south polar axis (if in the SPCH) in the EUVI FOV 
(i.e., events which occur at lower latitudes) show greater devia¬ 
tion in COR1 on average. This can be associated with the fact 
that the trajectory of the jets is not simply radial but bends to¬ 
wards the equator: the actual jet angle is the difference between 
the position vectors at the heights of 1 and 2 /? 0 (see the green 
dashed line that does not coincide with the red one). We can as¬ 
sume that jets propagate, on average, along the magnetic field 
lines, which are almost radial near the solar dipole axis, while 
those at lower latitudes deviate more markedly from the radial 
direction because of the dipolar structure. This is also consistent 
with the over-expansion towards low latitudes of the fast solar 
wind in polar coronal holes (e.g..lFis kl 1996b . The non-radial out¬ 


ward propagation of the jets is a property consistent with those of 
other coronal structure: e.g., coronal streamers and polar plumes 
extend non-radially (IWandll996l) . as well as non-radial seems 
to be the pro pagati on of CMEs, as found in earlier studies (e.g., 
ICremades & BothmerlUOOH) . It can be noticed that when going 
from the EUVI to COR1, the changes in PA are larger in the 
north coronal hole than in the south one. More precisely, linear 
fits of data points, with linfit. pro and linfitex.pro for the 
North and the South give the values of the parameter a and b that 
are summarized in Table [I] 

The parameter a, which represents the slope of the fitting line 
for the PAs, and hence a measure of the average bending of jets, 
is indicated as for the north (second column), and a s for the 
south (forth column) in Table Q] We also list the values of the 
parameter b, as b n for the North (third column) and b s for the 
South (fifth column), although they will not have a relevant role 
in the discussion. The bending is found to be on average greater 
in the North than in the South, giving an indication that jets 
are more deflected in the North than in the South {aulas > 1). 
The value of the quantity (a v - as)/as shows that is around 
26% larger than as . The square of the ratio of the coefficients, 
( 1 a^/as ) 2 , is also listed. The physical meaning of this quantity 
will be discussed in the last section. This asymmetry is also evi¬ 
dent in Fig. 12.II where we plot the absolute value of the relative 
jet-bending from the solar axis \0 COR1 - 9 EUVI \/\0 lt - Q EUV, \ as 
function of time for north polar jets (black void triangles) and 
south polar jets (red void squares). The parameter G„, which rep¬ 
resents the PA of the solar axis, is 0 deg at the North and 180 deg 
at the South. We can infer that there is not a particular temporal 
dependence of the PA displacement from EUVI to COR1 FOV 
on time, as might be the case for Ulysses measurements during 
the fast latitude scans, and the average magnetic deflection is 
larger at the North pole than at the South pole. 

The evidence that jets are more deflected in the north pole 
than in the south pole could be associated with the different 
bending of the open magnetic field lines: in our case, jets oc¬ 
curring in the North, are more deviated from their origin since 
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Table 1. Values of parameters for the linear fits of the EUV and COR1 PAs, latitudes, and longitudes for jets at the North and South poles. 



ay 

b N 

[deg] 

a s 

bs 

[deg] 

(ay - a s )/a s 

(a N /a s ) z 

PA 

LINFIT 

2.18+0.09 

-1.29 ±0.97 

1.72 ±0.09 

-129.20+ 15.44 

27% 

1.60 + 0.21 

LINFITEX 

2.25 ± 0.04 

-1.31 ±0.45 

1.78 ±0.04 

-138.93 + 6.98 

26% 

1.61+0.09 

Latitude 

LINFIT 

1.77+0.36 

-72.33 + 28.59 

1.45+0.24 

47.64 ± 18.26 

22% 

1.49 + 0.78 

LINFITEX 

2.31+0.06 

-115.51+4.97 

1.62 ±0.04 

56.74 ± 3.06 

43% 

2.04 + 0.15 

Longitude 

LINFIT 

1.11 ±0.07 

-10.06 + 4.66 

1.00 ±0.02 

-1.53 + 2.42 



LINFITEX 

1.04 ±0.01 

2.08+0.51 

0.98 ±0.01 

0.81 ±0.77 

- 

- 



Fig. 3. Absolute angular displacement normalized to the EUVI PA rel¬ 
ative to the solar axis as a function of time for jets occurred in the 
North (hollow black triangles) and in the South (hollow red squares). 
The quantity 6„ is 0 deg for events at the North and 180 deg for those 
at the South, in order to compare values between jets seen at the two 
poles. The dashed horizontal lines mark the average displacement: 2.34 
deg for jets at the North and 1.10 deg for those at the South. 


the magnetic field lines are more curved than in the South dur¬ 
ing this period of observations. 


2.2. 3D analysis 

The distinctive capability of STEREO is to perform 3D stereo¬ 
scopic geometry or t raje ctory reconstruction of solar structures, 
such as loops (e.g. Nistico et al.ll2013l). and CM Es (Bosman 
et al. 120121k Understanding the 3D evolution of jets through 
the corona as observed with the EUVI and COR1 instruments 
can provide additional information about their magnetic deflec¬ 
tion. In order to examine this aspect, we measured the 3D po¬ 
sition of jets in Stonyhurst longitude and latitude, at given ra¬ 
dial distances of 1 and 2 R Q . In practise, we used the routine 
scc_measure. pro, available from the SSW package, which al¬ 
lows to determine these quantities. We did this for a subset of the 
jets in the catalogue, i.e., for those events for which the determi¬ 
nation of the 3D position is more reliable because of better view¬ 
ing; this was possible for 38 events (20 at the North and 18 at the 
South, respectively). We identified the coordinates at 1 R e by tri¬ 
angulating the base of the jet in EUVI images. For each event we 
collected 10 measurements, in order to take into account errors 
due to the triangulation process, and calculated the average val¬ 
ues and standard deviations for the longitude and latitude. The 


same procedure has been performed for COR1, in order to mea¬ 
sure coordinates at 2 R Q . In this case, collecting points at this 
fixed distance required more efforts, since there is not any refer¬ 
ence that can help us in locating the jet (in the EUVI FOV, the 
footpoint of the jet or the limb is a good marker). The visibility 
of the jet, in both cases, has been eventually improved by using 
difference images. 

The results of the 3D measurement are shown in the top pan¬ 
els of Fig. [5] The jet location is de-projected in polar plots, show¬ 
ing the Stonyhurst longitude (concentric circles) and latitude (ra¬ 
dial lines) at 1 R Q as red dots, and at 2 R, :1 as green dots. A blue 
dashed line connects the locations of a jet, which is marked by a 
number according the catalogue from lNistico et al.l ( 2009 ). The 
radial lines and arcs centered on the dots are the error bars for 
the latitude and the longitude, respectively. It is worth noting that 
these graphs provide an anticipation of what we can see with 
Solar Orbiter, when it will be able to see directly the polar cap 
when orbiting out of the ecliptic plane. The first impression is 
the consistency of the measurements from EUVI and COR 1, 
which are taken independently, showing displacements toward 
low latitudes at higher distance from the Sun, and almost a ra¬ 
dial trajectory (i.e. a small longitudinal shift). In a similar way 
as done for the PAs, we can fit the measurements between 1 and 
2 R 0 for the latitudes and longitudes, respectively, by using 0. 
The scatter plots are shown in Fig. [5] The latitudes for both poles 
(top panels) are distributed far away from the bisector of the 
plane (dashed red line, which ideally should mark events with 
no change in latitudes), and the linear fits (in green with linfit 
and yellow with linfitex) return a slope a y = 1.78 - 2.31 for 
the North, and as = 1.44 - 1.62 for the South (see the second 
row group in Table [[]). In this case, the quantity (ay - as)/as 
shows that ay is around 22-43% larger than as, which is almost 
consistent with the results obtained from the PA measurements. 
The ratio (ay/a,s) 2 is around 1.5-2.0. On the contrary, the longi¬ 
tudes (third row group in Table Q} are very close to the bisector 
(slopes are 1.11-1.04 and 1.00-0.98 for the north and south jets, 
respectively), suggesting that the studied jets did not exhibit a 
significant shift in the azimuthal direction. 


3. A model for the coronal magnetic field 

From the results shown in the previous section, a natural ques¬ 
tion arises: “Is the different deflection of jets, in the North and in 
the South poles, an indication of a magnetic North-South asym¬ 
metry?”. There are several extrapolation methods for estimat¬ 
ing the magnetic configuration of the solar corona (for de tails 
see, I Altschule^&Newkjrrilj_969y KN As- 

chwanden 120051 chanter 5). In the potential field (i.e.. current- 
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Fig. 4. Polar plots of 3D positions of jets seen at the North (left) and South poles (right). The circles represent the heliographic longitudes, and the 
radial lines the latitude meridians measured in deg. Position of jets at 1 R 0 are in red, at 2 R 0 in green, respectively. The Sun-Earth direction is 
given by a red arrow in both plots. 


free) approximation, the magnetic field can be obtained as the 
gradient of a scalar potential, B = -V'F, and the scalar potential 
function can be expressed in spherical coordinates as an expan¬ 
sion in terms of the Legendre polynomials /""(cos 6): 

N l 

'fV, e, cp) = R e 2 2 f,(r)P ';\cos 6)[g? cos (m<p) + h'” sin (m<f>)]. 

1=1 m =o 

( 2 ) 

The indices Z and m are integer numbers: Z is strictly positive 
and defines the number of axes of symmetry of the field, while 
m can assume 2/ + 1 values (-Z,..., -1,0,1,..., Z) and defines the 
orientation of the axes in spherical geometry. The function: 


the dipole (Z = 1), the quadrupole (Z = 2), and the esapole (at 
1 = 3) terms, with the corresponding moments (or harmonic co¬ 
efficients) gi,g 2 , (D (we dropped the superscript m in the coeffi¬ 
cients since it is always 0). 

After some algebra, we find the magnetic field components 
due to the dipole: 



(4) 

(5) 


the components due to the quadrupole: 


Mr) = 


( r w /r) M 
0 r w /R B ) M 


( r/r w ) 1 
(Re/rJ’ 


(3) 


fixes the position of the solar wind source surface at r w solar 
radii; beyond r w the magnetic field lines are purely radial, re¬ 
producing the configuration of the magnetic field in the solar 
wind (iAschwandenll2005l) . Thus, Eqs.[2][3]must be used only for 
R Q < r < r w . 

The components of the magnetic field can be found as the 
derivative of the scalar potential <L. We can further simplify the 
expressions of the magnetic field components, if we assume ax¬ 
ial symmetry. This assumption is supported by the 3D analysis 
of the jet position, showing negligible shift in the longitudinal or 
azimuthal direction: then the component B :ll is null. Of course, 
this assumption implies that we are neglecting other effects like 
the possible magnetic dipole tilt. Axial symmetry is enforced by 
setting m - 0. Indeed, near the poles (6 - 90 deg, 180 deg), all 
associated Legendre polynomials P'” oc sin'” 6 —> 0, except for 
m = 0. 

The expansion now depends only on the index Z. If we trun¬ 
cate the series at Z = 3, we have three contributions that give 


B?\r,6) = I(^) 4 (^^!) 02 (3cos 2 0-1), (6) 

Bg\r,6) = 3^j ( ^'-^ 5 ) g2cosgsin6,; ( ? ) 

and the components due to the esapole: 

B?\r,6) = I(^| 5 (1^1L] 03 ( 5 cos 3 0-3cos0), (8) 

Bg\r,e ) = ^(—) (-y—03(15 cos 2 (ZsinO - 3 sin0)(9) 

^ ' r ' \ r w ~ RqI 

In the limit of r w —> oo, we have the classical expressions 
for the dipole, the quadrupole, and the esapole in free space. The 
resulting magnetic field can be written as the sum of the dipole, 
quadrupole, and esapole contributions: 

B r (r, 6) = B { r l) + B J. 2) + B?\ 

Bg(r,9 ) = B^+B^ + B^. 


( 10 ) 

( 11 ) 
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Fig. 5. Top: scatter plots of the latitudes mea¬ 
sured with EUVI at 1 R e , and COR1 at 2 R e for 
the north (left) and south (right) polar jets. Lin¬ 
ear fits of the points is given by the green and 
yellow dashed lines, while the red one is the 
bisector of the plane. Bottom: similarly to the 
previous graphs, scatter plots of the longitudes. 
The num bers in blue identify the events in the 
catalog of lNistico et all (2009). 


The magnetic field lines can be obtained by integrating, with 
a Runge-Kutta scheme of the 4th order for example, the two first 
order differential equations: 

I drjds — B r /B. 

\ ( 12 ) 

[ dd/ds = Bg/rB. 

where s is the distance along the field line. This model is then 
used to find a relationship between the coefficients a and the nor¬ 
malised multipole coefficients of the reduced field model. 

3.1. Magnetic moments from the Wiicox Solar Observatory 

Given the model for the coronal magnetic field with the as¬ 
sumption of axial symmetry and the truncation to the dipole, 
quadrupole and esapole moments, we can search which val¬ 
ues of the magnetic moments gi, go, gs, are suitable for better 
describing the jet magnetic deflections. For comparison, these 
coefficients are calculated from magnetograms data, provided 
and published by the Wilcox Solar Observatory (WSO) (see the 
website http: //wso. stanford.edu/). For inferring values of 
these coefficients from magnetograms, two kind of hypotheses 
are made regarding the inner bo undary conditions in the pho¬ 
tosphere (IWang & Sheelevlll992 ). Indeed, from magnetograms 
we can measure the line-of-sight (LOS) component of the pho- 
tospheric field at a given latitude a on the solar disk. The “clas¬ 
sic” model takes into account the projection of the LOS pho- 
tospheric field along the radial and latitudinal component ( B, = 
Bros cos a , Bg = Bios sin a , with a the line-of-sight angle); con¬ 
versely, the “radial" model assumes that the photospheric field 


is totally radial in the photosphere ( B , = Bios / cos a, Bg = 0). 
Since the ma gnetic field is non poten tial and nearly radial at the 
photosphere (IWang & Sheelevl[l992 ). a better approach is con¬ 
sidered to be the radial model. 

The top panels of Fig. [ 6 ] show the temporal evolution of the 
coefficients for the “classic” and the “radial” model, as calcu¬ 
lated by the Wilcox Solar Observatory. The bottom plots are the 
temporal evolution of the ratios g/g™ and gY/gY (where the 
superscript W signifies g coefficients calculated by the WSO). 
The region bounded by the dashed lines is the temporal window 
in which our jets are observed. 

We can see that the harmonic coefficients evolve in time, 
and that the quadrupole moment is less in magnitude than the 
other ones and exhibits an oscillatory behaviour, especially in 
the radial model, with alternating sign and with a period of ap¬ 
proximately two year s, indicating a possible association with the 
biennal oscillations ( Vecc hio & Carbonel 120081 : IVecchio et al.1 
l2012l : lBazilevskava et alJl2014l) . On the contrary, the dipole and 
esapole moments show a constant sign. Below, we try to calcu¬ 
late these coefficients in a different way and to compare them 
with those obtained from the WSO. Further, we find the asso¬ 
ciated coronal magnetic field structure and the position of the 
heliospheric current sheet projected on the solar surface. 


3.2. Fitting the magnetic field model to the polar jets PAs 

Here, we estimate the best values of the coefficients <72 and <73 
that fit our jet observations and measurements of PA, in terms 
of the dipole moment (which is used as a normalization factor). 
For several values of §2 = ( 72/(71 and §3 = ( 73/(71 coefficients, 
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Harmonic coefficients (classic model) Harmonic coefficients (radial model) 




From 2003-11 -19 to 2010-08-09 

Fig. 6. Top: Time evolution of the harmonic coefficients g^, g*,g^, as calculated by the Wilcox Solar Observatory (http: / /wso. Stanford. edu) 
in the classic (left) and radial (right) model. Bottom: time evolution of the quadrupole and esapole moments normalised to the dipole component. 
Note the periodic change of the quadrupole over a scale of 20 Carrington rotation (~2 yr) in the radial model. The vertical dashed lines enclose 
the period of our observations. 


we integrate numerically the equations (IT2l) from the base of the 
jets, z'.e.the PA as measured in the EUVI FOV, until 2 R Q . This 
yields the difference between the final PA from the numerical 
integration, g CORl Mo^ w hich is a function of <72 and fa, and that 
one measured in the COR1 FOV, (f OR1 . Then we calculate the 

standard deviation cr(fa, § 3 ) as: 


fafa, fa) = \ 


n 


i=i 


ef 0 R 1 ^(fa,fa)- 


■ 6 


t CORlobs 


N- 1 


(13) 


with N the total number of magnetic field lines successfully in¬ 
tegrated from the jet base at 1 R Q up to 2 R Q . The number N 
is not necessarily 79 but can be less, depending on the values 
of the magnetic moments since the magnetic field lines can be 
closed without reaching 2 R Q . Fig. [7] shows for example some 
magnetic field lines integrated for a few jets with given values of 
the coefficients fa and fa. We can notice that in some events a 
good agreement between the final position from the integration 
and the observed PA at 2 R Q is found, some others show a con¬ 
siderable gap, whilst for an event the integrated line results to be 
closed and does not reach 2 R Q . The coefficient fa spans from -3 
to 3 and fa from 0 to 3, both in steps of 0.1. We give the results 
of cr(fa,fa) as 2D contour maps representing the value of the 
standard deviation as a function of the quadrupole (vertical axis) 
and esapole (horizontal axis) moments, using different models of 
the magnetic field, i.e ., classic and radial, and different distances 
of the solar source surface r,„ = 2.5 and 3.25 R.. (Altschuler & 
Newkirk ! 1969fc iHoeksemall 19951) . The grid in the maps is deter¬ 
mined by the varying values of fa and fa, with a resolution of 
0 . 1 for both parameters, as used in the numerical model. 

The top panels of Fig. [ 8 ] show the standard deviation maps 
for the classic (left) and radial (right) model considering the so¬ 


lar source surface at 2.5 R Q \ the bottom row gives the same re¬ 
sults for r w - 3.25 R 0 . The colorbar provides the variation range 
of cr: from 4 (black) to 16 deg (red). The purpose of these maps 
is to give us an indication of which values of fa and fa minimize 
the standard deviation (dark regions in the maps), and better fit 
the deflections of jets. In addition, they are overlaid by some 
coloured lines that enclose portions of the maps characterised 
by the same percentage of jets, hence magnetic field lines, suc¬ 
cessfully integrated. Indeed, if a certain value of cr is obtained 
for different values of the magnetic moments, the number of in¬ 
tegrated lines N will provide a further constrain to discriminate 
which values of fa and fa are more reliable: the higher N, the 
better the adaptation of the magnetic field model to our jet ob¬ 
servations. For this reason, we can exclude the minimum in cr 

found for fa -2.5 and fa ~ 1.8, since it is obtained for as few 

as almost 50% of the observed jets. 

We do not obtain specific and exclusive values of magnetic 
moments that fit our observations. The shape of the low-cr is 
almost similar in all cases and appears to be more sensitive to 
the esapole moment, since it is more narrowed for some values 
of fa, and spread for several values of fa, but a patch with cr < 6 
deg and in the limit of N — 90% is found for negative values of 
the quadrupole, in agreement with the values from WSO. This is 

marked by red crosses in the maps at values of (fa) -0.5 and 

( fa )~ 1.6 — 2 . 0 . 

We can now infer the structure of the coronal magnetic field 
by computing Eqs. (4-9) with the obtained values of (fa) and 
(fa) and plotting the magnetic field lines. This is shown in Fig. 
[9] the HCS in the considered period of observations results to 
be coned southward, forming an angle of about 10 deg, which 
is broadly consistent, although somewhat larger, with some es¬ 
timates found in the literature ranging between 3 and 10 deg 
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Fig. 7. Example of magnetic field line integration for some jets with 
moments g 2 = -0.8 and g 2 = 1.3. The starting points of the integration 
are at 1 R e from EUVI observations and represented as coloured dots. 
The integrations is made up to 2 R Q and the results (void squares ) are 
compared with the corresponding jet positions as observed with COR1 
(plus signs). For some events, especially those at lower latitudes, the 
integration can result in a closed line (purple dashed line). 


jSjmjrson^talJl996^ Heber et al.ll9%l : lMursula & Hiltulal2004l : 


lErdos & Baloghl 20101) . On the other hand, the HCS is not a sta¬ 

tionary feature and the tilt is subject to a change over the time 
due to the evolution of the magnetic field structure. For the pe¬ 
riod under examination, the average values of the quadrupole and 
esapole components, as determined by the Wilcox Solar Obser¬ 
vatory, are of the order of g^/g^ ~ -0.1,0.05 and g™/g™ ~ 
1.0,1.6 for the classic and radial model, respectively. In addition 
the quadrupole component shows a large variability, assuming 
positive and negative values with peaks at ±0.5c/j l in the case of 
the radial approximation. 


4. Discussion and conclusions 


In this work, we used polar corona jets as a probe for understand¬ 
ing the magnetic field structure of the solar corona during a solar 
minimum. Since the corona is an environment with a plasma-/! 
parameter less than 1, we assumed that jets on average propagate 
along the magnetic field lines. For simplicity, we also assumed 
that the large scale solar magnetic field during a solar minimum 
is axisymmetric. We measured the PA of jets at 1 R 0 in the EUVI 
FOV, and at 2 R 0 in the COR1 FOV, and analysed the deflection 
of jets. We found that jets are deflected toward low latitudes, in 
agreement with the fact that the magnetic field lines are bent to¬ 
ward the equator, and this deflection is more pronounced in the 
North pole than in the South pole. This North-South asymme¬ 
try has been found in other d atasets, start ing from photospheric 
magnetic field measurements (lHo eksemall l995), glo bal distribu¬ 
tion of the solar wind speed (iTokumaru et al.il2015l). anal ysis o f 
latitudinal gradient in s olar energetic particles dsimpson et al.l 
ll996klHeber et al.ll 19961) . and also magnetic field measurements 
in the interplanetary medium by the Ulysses spacecraft (Erdos 



Fig. 8. Contour maps representing the_standard deviation (cr) from ( 1131 
for different values of r w (top: 2.5 R 0 , bottom: 3.25 R a ), for the classic 
(left) and ndial (right) model, as a function of g 2 (vertical axis) and g 2 
axis) coefficients. The colorbar defines the variation rang; of 
on, coloured lines enclose regions with different percentage 
(100, 95, 9), ...) of jets, hence magnetic field lines, which are succtss- 
fully integrated. The red crosses locate the chosen minimum in a low-cr 
at least 90% of jets integrated. This minimum is found for 
. and (g 2 ) = 1.6 (bottom panels)- 2.0 (top panels). 


(horizontal 
cr. In additi 


region wit! 

(g 2 ) = -o. 


& Balogh l20ld) . This asymmetry can be modelled in terms of 
multipole compo nents of the global magnetic held ( Bravo & 
Gonzalez-Esparza lSoOOl i lMursula & Hiltulall2004l) : during solar 
minima the quadrupole moment tends to influence the total mag¬ 
netic field, weakening it in one pole and enhancing it in the oppo¬ 
site one. We have estimated the contribution of the quadrupole 
moment, starting from jet PA data, and comparing our results 
with those of the WSO. We expressed the coronal magnetic field 
as the sum of the dipole, quadrupole, and esapole moments, start¬ 
ing from the expression of a scalar potential function <t> in terms 
of spherical harmonic expansion. We traced magnetic field lines 
from the footpoint of jets (at 1 R Q ) until 2 R Q for different values 
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Fig. 9. Structure of the coronal magnetic field lines with (g 2 ) = -0.5 
and (§ 3 > = 1.6, and the corresponding HCS position projected on the 
source surface at 3.25 R e for the interest period of observations between 
March 2007 and April 2008. Inward magnetic field lines are plotted in 
red, outward in blu. 


of the magnetic moments. We obtained the harmonic coefficients 
g 2 and § 3 , normalised to the dipole, that minimize the standard 
deviation of position angle differences coming from the numeri¬ 
cal simulation and the observations. We obtained as reliable val¬ 
ues (g 2 ) = -0.5 and (g 2 ) = 1.6. From these estimates, we were 
able to compute the structure of the magnetic field lines (Fig. [9]) 
in which the heliospheric current sheet (HCS) shows an offset 
of about 10 deg, a value c onsistent with the results obtained by 
iMursula & Hiltul al (120041). but somewhat larger than that found 
bv lErdos & Baloghl ( 201 Oh . whose observations however corre¬ 
spond in part to different periods and are taken at much larger 
distances from the Sun, around 1.4 AU during the fast latitude 
scans of Ulysses. 


In our analysis, we show that the slopes of the linear fits for 
the PA and the latitudes (from 3D measurements) are different 
between the two poles. The angular coefficients of the fits can 
be immediately related to the ratio of the magnetic field val¬ 
ues ZN IB m — 1.12 - 1 . 21 , as reported in Table 1 from Erdos 
& Balogh (1201 Oh . according to the following geometric inter¬ 
pretation in terms of conservation of magnetic fluxes through 
the Sun’s poles. Consider the sketch of the Sun in Fig. [TOj with 
the polar caps marked by dashed lines at the distance of 1 R Q , 
which can essentially represent areas embedded in open field re¬ 
gions, such as the polar coronal holes. We can express the mag¬ 
netic fluxes <t> approximately as the product of the area A of 
the polar cap with an average polar magnetic field B. The sur¬ 
face of a spherical cap depends on the half opening angle 6 as 
A - 2nR^{\ - cos 6 ), which we assume equal for both poles at 1 
R e , for example. The angle 9 is analogous to the PA in our mea¬ 
surements. Hence, because of the asymmetry, the flux Os in the 
South pole will be larger than O v in the North. Thus, we have at 


the distance of 1 R a : 


An(Rq) 


A s (R q ) 


Oy(Ao) 

O s (fl 0 ) 


Rjv(Ro) 

Bs(Re) 


(14) 


Moving away of the Sun, the magnetic field diminishes with 
the distance r and the area of the projected polar cap must in¬ 
crease because of the flux conservation ( 0 ,(R G ) = 0 ,(r) with 
i = N,S). We can consider the variation of the magnetic field 
normalised to the distance Z? G as B,(r) = B i (r)B,{Rr,A and the 
expansion of the polar caps with the distance can be addressed 
in terms of a variation of the opening angle 6 of a factor a, which 
is different between the two hemispheres (and also depends on 
the radial distance). Thus, at a certain distance r (which is taken 
to be R Q < r < r w in the PFSS model) the polar cap areas can be 
expressed as: 

Aj(r) = 2nr 2 [1 - cos (a,(r)0)] (15) 


with; = N, S. A second order approximation for the cosine func¬ 
tion gives cos (a9) w 1 - \a 2 0 2 , and, finally, the area can be ex¬ 
pressed as A,-(r) = nr 2 (a 2 9 2 ). By taking into account the ratio of 
the magnetic fluxes between North and South at r, we have: 

<&N(r) _ B N (r)A N (r ) ^ B N (r) B N (R 0 ) / a N (r) \ 2 
®sM B s (r)A s (r ) ~ B s (r) B S (R Q ) \ a s (r) j 

and by combining eqq. (IT4l) and (IT6l) in virtue of the mag¬ 
netic flux conservation, we obtain the final relation: 


Bs(r) ~ /'flvW'T 
B N (r) UsO)/ 


(17) 


which links the estimated deflections with the ratio of the mag¬ 
netic field magnitudes. The squared values for the ratio a^/as 
at r — 2 Rq, according our analysis of the jet deviations, are 
around 1.5-2.0 (see Table[l]i, almost 25-65% larger than the ra¬ 
tio of the magnetic field estimated bv lErdos & Baloghl (l2010i) . in 
agreement with the larger estimate of the HCS offset in our anal¬ 
ysis. Theagreementbetween^ourresult^andthose^rom Erdos & 
Balogh T 2010 l) are rather satisfactory. We would like to outline 
the comparison is made on the basis of different observable (jet 
deflection against interplanetary magnetic field measurements 
normalised to 1 AU) and refer in part to different periods on time. 
In addition, a list from the Wilcox Solar Observatory of the tilt of 
the HCS reports values greater than 10 deg for the period under 
interest. 

Therefore, we have an independent indication that the solar 
magnetic field can indeed exhibit a N-S asymmetry, a result that 
can have profound implications on the models of solar dynamo. 
The forthcoming missions Solar Probe Plus and Solar Orbiter 
will have a crucial role in assembling more precisely the extent 
of the N-S asymmetry, thanks to both in situ measurements of the 
nearly coronal magnetic field and remote observations. In partic¬ 
ular, the Solar Orbiter UV instruments will allow to accurately 
check the size of the polar coronal holes and the deflection of 
polar jets, allowing to better constrain the solar magnetic field. 
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Fig. 10. Sketch representing the N-S asymmetry of the solar magnetic 
field in terms of magnetic fluxes through the polar caps. The polar caps 
on the solar disk at the distance of 1 R a represent ideally regions in polar 
coronal holes with equal area, which depend upon the half-opening an¬ 
gle 9. Because of the asymmetric over-expansion of the magnetic field 
lines between the two hemispheres (represented by the curved dashed 
lines departing from the solar disk), the angle 9 at a given distance r 
will increase of a factor an in the North and as in the South, with 
fljv > Os- The resulting size of the Northern polar cap will be larger 
than the Southern one. 
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